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We address the question to what extent JIMWLK evolution is capable of taking into account 
angular correlations in a high energy hadronic wave function. Our conclusion is that angular (and 
indeed other) correlations in the wave function cannot be reliably calculated without taking into ac- 
count Pomeron loops in the evolution. As an example we study numerically the energy evolution of 
angular correlations between dipole scattering amplitudes in the framework of the large Nc approx- 
imation to JIMWLK evolution (the '"projectile dipole model"'). Target correlations are introduced 
via averaging over (isotropic) ensemble of anisotropic initial conditions. We find that correlations 
disappear very quickly with rapidity even inside the saturation radius. This is in accordance with 
our physical picture of JIMWLK evolution. The actual correlations inside the saturation radius in 
the target QCD wave function, on the other hand should remain sizable at any rapidity. 



I. INTRODUCTION AND CONCLUSIONS 

The CMS observation of angular and long range rapidity correlations in the hadron sgectrum^he so called "ridge" 
in proton-proton collisions [1], has triggered a lot of discussions in recent literature [5] , [3] , 3 , [H| • In particular the 
approaches of the three last papers [3],[^],[5l] are based on the idea that correlated gluon emission is due to the 
correlations in impact parameter plane preexisting in the incoming wave function of the target and projectile hadrons. 
Such correlations certainly exist in a hadron wave function, and in the context of high energy evolution can be encoded 
in the initial conditions for the evolution. It was also argued in [4] that the correlations are leading effect in The 
purpose of the present note is to address the question what is the fate of such correlations as the hadron is evolved 
to high energy. In particular we ask whether these correlations can be studied by evolving the target/projectile wave 
functions with JIMWLK evolution [6] (or the Balitsky-Kovchegov (BK) equation, which is its large Nc limit [7]). To 
this end we perform simple numerical calculations and supplement them with qualitative analysis based on physics of 
the JIMWLK evolution. Our numerics is performed in the framework of the dipole model approximation to JIMWLK 
evolution, and thus is sensitive only to the leading Nc part of JIMWLK, but we believe that, with minor modifications 
our conclusions are valid for full JIMWLK evolution as well. 

Our conclusions are the following. We find that the JIMWLK evolution leads to exponentially quick disappearance 
of correlations (including angular correlations relevant for gluon emission) with rapidity. This disappearance is 
straightforward to understand. It is the consequence of the fact that JIMWLK evolution is valid only for color field 
modes with transverse momenta smaller than the saturation momentum where one indeed does not expect correlations 
to be present. As was discussed in [4], we expect correlations to be present only for points in transverse plane within 
the saturation radius of each other, and therefore for momentum modes greater or equal to the saturation momentum. 
The evolution of these modes even in a dense hadronic wave function is not governed by the JIMWLK evolution but 
rather by KLWMIJ evolution [9] for k ^ Qs and by the Reggeon Field Theory (RFT) including Pomeron loops jlQ| 
for k ^ Qs- 

The failure of JIMWLK to properly account for correlations should be understood in the following way. The evolved 
wave function of a hadron does indeed contain correlations in impact parameter space, even in the leading order in 
expansion. This has been shown analytically and numerically for angular independent correlations within Mueller's 
dipole model in [l2|, and there is every reason to expect that this is also the case for angle dependent correlations 
of interest to us. However JIMWLK evolution approximates the scattering amplitude of two dipoles on a hadronic 
target by contributions where the two dipoles scatter on gluons with vastly different rapidities (this has been dubbed 
"'long range multiple scatterings"' in [13]). Since the correlations in the wave functions are between gluons which are 
close in rapidity to each other, these correlations are simply not included in the JIMWLK equation. While JIMWLK 
approximation does properly account for leading contributions to the scattering amplitude in the saturated regime, it 
significantly underestimates those for small (smaller than saturation radius) dipoles. To account properly for rapidity 
evolution of the scattering amplitude of small dipoles, one needs to evolve it with KLWMIJ evolution equation. The 
KLWMIJ evolution does indeed include the "'short range multiple scattering contribution" [isj and should correctly 
describe correlations in impact parameter plane. The obvious complication here is, that since we are interested in 
a dense target (corresponding to high multiplicity events in the CMS data), it certainly contains low momentum 
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modes which evolve according to JIMWLK. Thus we are faced with the situation where proper treatment of angular 
correlations at high energy requires us to include both KLWMIJ and JIMWLK evolutions within the same framework, 
and in this sense we have to deal explicitly with the Pomeron loop effects. In fact this necessity is even more acute, 
since the Pomeron loops give leading contributions to the evolution of the modes at A: Qs, and it is presumably 
these modes that contribute most to angular correlations. 

We thus conclude that future attempts to properly numerically estimate the size of correlations at high energy will 
require explicit inclusion of the Pomeron loop effects. We note that the notion that it is the Pomeron loops that 
are crucial for correlations at high energy is not new [14], here we merely recast this argument in the framework of 
JIMWLK/KLWMIJ evolution. 

The structure of this paper is the following. In Sec. 2 we recap the arguments of [4] about angular correlations in 
gluon emission, recasting them in a somewhat more transparent semiclassical form. In Sec. 3 we present results of 
our numerical calculations of the evolution of angular correlations in the dipole model (leading Nc JIMWLK). Finally 
in Sec. 4 we discuss the interpretation of these results based on the physical picture of JIMWLK/KLWMIJ evolution 
and flesh out the arguments for necessity of Pomeron loops. 



II. ANGULAR CORRELATIONS 

In our previous note we discussed a simple picture of long range rapidity correlations and angular correlations 
between particles produced in a collision of two high energy dense objects. This quialitative picture also underlies 
the calculations of , @ . Long range rapidity correlation is an almost trivial consequence of boost invariance of a 
projectile wave function at high energy. Consider high energy scattering of a hadronic projectile on a stationary target 
in the lab frame. In the lab frame, the incoming particles are very energetic and they scatter by a very small angle 
with p+ ^ Pt- Thus recoil is negligible and eikonal approximation is applicable at high enough energy. Since the 
projectile is very energetic, its wave function is approximately boost invariant. In a boost invariant wave function 
gluon distribution at rapidity Yi and Y2 are the same. These gluons scatter on exactly the same target, and thus 
whatever happens at Yi also happens at Y2. If for a particular target field configuration a gluon is likely to be produced 
at Yi at some impact parameter, a gluon is also likely to be produced at Y2 at the same impact parameter, thus leading 
to long range rapidity correlations. Thus the long range rapidity correlations come practically for free whenever the 
energy is high enough so that the wave function of the incoming hadron is approximately boost invariant, and there 
is very little in the actual dynamics of the collision that can affect this feature. 

To understand why angular correlations also naturally arise in the context of high energy let us briefly recap our 
understanding of the transverse structure of the hadron in the saturation regime. It is convenient to think of the 
distribution of the (color) electric field configurations in the target. 

The target wave function is characterized by the saturation momentum Qg. The saturation momentum plays 
a dual role in the hadronic wave function. First, it measures the typical magnitude of electric field in the wave 
function. The scattering amplitude of a dipole on the target is given in terms of simple parton scattering amplitude 

S{x) = j as N{r) = 1 - -^Tr[S\{))S{r)]. The vector potential is simply related to the electric field 

as diA~ = F~\ Let us for convenience define electric field integrated over the longitudinal extent of the target, 
Ei = J dx^F~^ . The dipole scattering amplitude is then given in terms of gE^ and assuming for illustrative purposes 
that odd powers of E average to zero in the hadronic ensemble, we have roughly 

Ar(r)-l-e-4(^^"'-^)' (2.1) 

This is of order unity for = Q~'^ = {gE)~'^. 

On the other hand it is known that the field components with transverse momenta pt < Qs are suppressed in 
the wave function [15]. This means that the electric fields in the target are correlated on the length scale A ~ Q~^^ • 
Thus the saturation momentum doubles up as the inverse of the correlation length of target color fields. Typical field 
configurations in the target can thus be thought of having a domain like structure of Fig. 1. 

Now consider a projectile parton with charge q impinging on one of the domains of the target. While traversing 
the target field, the parton acquires transverse momentum 

^P = gq J dx^F- = gqE. (2.2) 

A parton at a different rapidity but with the same charge will pick up exactly the same transverse momentum if it 
scatters on the same "domain". This of course results in positive angular correlation of produced gluons. 

We note in passing that this simple picture also explains the fact noted in [J that angular correlations at angle (j) 
and (j)-\-7r have equal strength. At high energy particle production is dominated by gluons. Gluons of course belong to 
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FIG. 1: Typical color electric field configuration in the target. 

real representation of the gauge group, thus it is equally probable to find an incoming gluon with charge q and charge 
—q in the projectile wave function at any rapidity. Suppose, for example that a given configuration the color field in 
the target is in the third direction in color space = EiS^^^ while the gluons in the incoming projectile corresponds 
to the vector potential in the second direction Af. One can always write = — z/2(A — A*), where A = Ai -\- iA2 is 
positively charged with respect to color charge in the third direction, and A~ is negatively charged. Thus necessarily 
equal number of gluons in the incoming projectile have opposite sign charges and are kicked in opposite directions 
while scattering on the target. This produces equal strength correlations at angles zero and tt. This is not the case 
for quarks which carry fundamental charges, and it is quite clear that taking into account the projectile quarks will 
lead to stronger positive angular correlation than the negative one. 

Going beyond the qualitative picture described above, the two gluon inclusive production probability discussed in 
(3 is given by [ill (see also [13) [Si 



dPpdPkdrjd^^ 
with 



dN 

= {a{k)a{p))p,T (2.3) 



cT{k) = I e^^^^-^") / f{-z - xi) • /(xi - z) P{xr)[S\x^) - S\z)][S{x^) - S{z)] p{x,) . (2.4) 



Here 



= fcf (2.5) 

and p = —iT^p^. In these formulae p^{x) is the valence color charge density in the projectile wave function, while 
S^^{x) is the eikonal scattering matrix determined by the target color fields. The charge density is normalized such 
that for a single gluon p^ = gT^. The two gluons here are produced independently of each other, but from exactly 
the same configuration of color charge sources while scattering on the same target field configuration. 

The average in eq. ()2.3p denotes averaging over the projectile and the target wave functions. The averaging over p 
is understood as averaging over a classical ensemble with a probability distribution function Wp [p] [6l] 

{0)p = j DpWp[p]0 . (2.6) 

and similarly for the target average. After averaging over all target and projectile configurations, the single gluon 
emission amplitude (<j) p^t must be isotropic. However for any given configuration it is not isotropic and peaked in one 
particular direction. This anisotropy produces angular correlation among the emitted gluons as discussed above. To 
reiterate, for a fixed configuration of the projectile sources p{x) and target fields S{x)^ the function (j{k) as a function 
of momentum has a maximum at some value k = q. Therefore the product in eq. (|2.4p is maximal for k = p = q. 
The value of the vector q of course differs from one configuration to another, but the fact that momenta k and p 
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are parallel does not . Therefore after averaging over the ensemble dPN/dkdp has maximum at relative zero angle 
between the two momenta (as we have explained above, there is actually the second degenerate maximum at relative 
angle = tt). The strength of the maximum of course depends on the detailed nature of the field configurations 
constituting the two ensembles (the projectile and the target). 

Thus, angular correlations emerge as a result of target/projectile averaging procedure over isotropic ensembles of 
anisotropic configurations. 

We note that e£. (|2.4p holds in the case when one of the colliding objects is dense and another one is 
dilute [131, [l85,[l9]. This is most likely not quite the situation encountered in the high multiplicity p-p events 
at LHC, where the density in the proton wave function is likely still not parametrically large, but is already not 
perturbatively small. The main features of our discussion are however born out by the expression eq. (|2.3p and we 
believe this approximation to be qualitatively correct. We note that the numerical calculations of [3^] use a further 
perturbatively expanded version of eq. (|2.4p . The approach to the dense-dense regime has been developed in [20| and 
has been used in [21]. It is however not clear that in the LHC environment it is quantitatively more reliable than 
simple expression eq. (|2.4p . When the dense system is produced in the final state, we expect the correlations produced 
by this mechanism to be washed out by the final state interactions and finally disappear for a very dense final state. 
There may be an alternative mechanism of producing angular correlations via radial fiow effects [22] which could be 
relevant to the ridge structure observed at RHIC[23], but this is far beyond the scope of thepresent work. 

Returning to eq. (|2.3p . we note that angular correlations should be the leading 1/Nc effect [J IH, [111 . The leading 
Nc piece in eq. ()2.3p comes from the configuration where the charge densities in each one of the single gluon production 
amplitudes are in the color singlet. The relevant average to calculate is 

{p''{xi)p''(xi)p\x2)p\x2))p{Ty{[S\xi) - S\z)][S{xi) - S{z)\}Ty {[S\x2) - S\u)\[S{x2) - S{u)\})t ■ (2.7) 

On the target side one needs to calculate averages of observables of the type described in the large Nc limit by the 
dipole model ^] 

{TT{[S\x)S{z)]}TT{[SHy)S{u)]})T = {s{x,z)s{z,x)s{y,u)siu,y))T (2.8) 

where s{x^y) = -^TT[Sp{x)SF{y)] - is the scattering amplitude of the fundamental dipole, and the equality in eq. (j2.8p 
holds in the large Nc limit. The approximation which is frequently used in the literature to calculate the averages of 
this type also invokes factorization 

{s{x, y) s{u, v)) = {s{x, y)){s{u, v)) . (2.9) 

Strict factorization of the type eq. ()2.9p is only possible if the statistical ensemble either consists of a single configu- 
ration, or the target fields on which the two dipoles scatter are completely independent of each other. There is of 
course no reason to expect that in the large Nc limit fiuctuations around some leading configurations are suppressed 
by powers of 1/Nc. Likewise, since we are interested in dipoles which scatter within correlation radius of each other, 
the field configurations should be by definition correlated. Thus the factorization eq. (j2.9p is not appropriate for study 
of correlations. 

Our objective in this note is a pilot study of the evolution in energy of correlations which one can encode in the 
initial target ensemble. As a framework for the evolution we take the projectile dipole model, which describes the 
evolution of s in the leading Nc limit. It is important at this point to avoid confusion and understand clearly which 
dipole model we are talking about, as there are two distinct approximations to the high energy evolution which are 
both sometimes called dipole model. 

In the first, which we will call target dipole model, one follows the evolution of the target wave function in terms 
of the density of dipoles (and their cumulants). This evolution as formulated in[ll| does not take into account finite 
density effects in the target wave function and is the large Nc limit of the KLWMIJ evolution[9]. Ref.[12] studied 
correlations and fiuctuations of dipole density in this approach for a single dipole target and found them to be 
significant in the leading order in 1/Nc. 

Another dipole approach, which we will refer to as the projectile dipole model, evolves the projectile wave function 
according to dipole evolution. The projectile scattering amplitude is then calculated by approximating the scattering 
amplitude of each projectile dipole by an eikonal factor. This approximation can be reformulated as the evolution of 
the target wave function. In this form it is a large Nc approximation to the JIMWLK evolution of the target wave 
function, that is, it does indeed take into account nonlinearities in the target evolution. 

The two dipole approximations implement very different physics in the target wave function. Our choice of the 
projectile dipole approximation is motivated by the expectation that high density effects in the target evolution should 
be important, and also by the fact that it is this approximation (or JIMWLK which includes corrections to it 

[26| ) that is used in current numerical studies of ridge [3]. We note however that this approximation does not take into 
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account proper splittings of the target dipoles, as stressed in [27]. As we show in the next section numericahy, and 
explain qualitatively in Sec. 4 this deficiency turns out to be crucial in the inability of this approximation to correctly 
evolve correlations present in the initial ensemble. 

The probability distribution of the dipole model W[s] evolves with rapidity according to [sl. [25j 

with a - the 'tHooft coupling, which is finite at infinite Nc. Our strategy is to choose an ensemble Wq[s\ of initial 
configurations s{x^y)^ which contains nontrivial angular correlations. Each configuration of the ensemble is evolved 
independently according to the BK equation [8] . The correlations at the final rapidity are then calculated by averaging 
the correlator over the ensemble of solutions 5y(x,?/). 

j DsWY[s\s{x,y)s(u,v) = j DsWQ[s\sY{x,y)sY(u,v) (2.11) 

where SY{x^y) is the solution of the BK equation with initial condition s{x,y). 

This procedure is similar to the one implemented in [28]. However, the focus of [l^ was in fiuctuations of the 
saturation scale and thus all configurations in the initial ensemble in [23] were chosen to be isotropic. In order to 
study angular correlations we have to allow the individual members s{x — y) of the initial ensemble to be anisotropic. 
The rotational invariance is restored by averaging over the whole ensemble rather than configuration by configuration. 

Here we report on our initial results, which mostly aim at qualitative understanding of rapidity dependence of 
angular correlations within the projectile dipole evolution. Like in most BK studies we do not consider impact 
parameter dependence. We view the resulting correlations as correlations at a fixed impact parameter, thus strictly 
speaking study of correlations in impact parameter plane are beyond our current calculation. Nevertheless we do 
not expect our results on weakening of correlations with rapidity at a fixed impact parameter to be affected by 
configuration by configuration fiuctuations in the impact parameter plane. 

For the sake of simplicity we do not calculate the two gluon production rate eq. ()2.3p but rather examine the simplest 
observable that can exhibit angular correlations - the correlator of two dipole scattering amplitudes s{x^y)s{u^v). 



III. BK EQUATION, INITIAL CONDITIONS, AND ANGULAR DEPENDENCE 

A. BK equation and initial conditions 

The BK equation for the imaginary part of the dipole scattering amplitude N{f) = 1 — s{f) (assuming impact 
parameter independent configurations) is: 

dYN{r) = %^ / ^V— [N{f) + N{r-f) - N{r) - N{r') N{r - f")] . (3.1) 

Z TT J f ^ (^f — fy 

Here f=x — yis a vector (in the transverse plane) connecting the two legs of the dipole, and r = |r|. 

As discussed above, we have to specify the initial ensemble of configurations of N{f). We choose all the configurations 
of the ensemble to have similar radial dependence, but distribute them homogeneously with respect to angle 6. A 
representative configuration at some initial rapidity Yq = \nl/xo is taken as 

7V(ro,rl = 1 - Expl^ - ar^xg^o<'^^<^Hxo,4/r')Fi0)^; a = . (3.2) 

Apart from the angle dependent function F, this is the same initial condition as used in f29| to fit HERA data. In 
1 29], the F2 low X data were reproduced using the BK equation with running coupling [32[. The parameters used in 
[29|] are: xq = 10~^ and the effective proton's radius R fitted to the F2 data, = 3.1 {GeV~'^). 

The function F{6) takes into account angular modulations of the scattering amplitude relative to some axis thereby 
refiecting anisotropy of a given target field configuration. For our study we choose F{6) of the form: 

m = 1 + lcos\0). (3.3) 

In contrast to previous numerous studies of the BK equation, in which the initial conditions were parameterized with 
respect to the dipole size only, eq. (j3.2p provides 2-dimensional initial data set (Fig. [2l For all the plots r is given in 
GeV-i). 
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FIG. 2: Initial conditions (|3.2|) . Left: profile of A/" as a function of r and 0. Right: the same in polar coordinates. 




L.S L.S L.S 



FIG. 3: as a function of r and at various values of rapidity: Y — Yq ^4.6, y = 6, Y = 10 

The cos^iO) dependence is motivated by dipole interaction with a constant target chromo-electric field eq. ()2.ip . 
In fact, even though this is quantitatively not quite true, qualitatively the scattering amplitude (|3.2p can be thought 
of as the scattering amplitude on a fixed, constant chromo-electric field configuration in the target. The function F{0) 
has period tt and is symmetric under 9 tt — 0. For this reason in what follows we will be quoting only results in 
the first quarter. 

We now define the ensemble of initial conditions by homogeneously distributing the direction of the target field in 
the impact parameter plane. In practical terms this amounts to shifting the angle in eq. (j3.3p by another angle 5, 
which is taken to be a random variable with the probability distribution W[6] = l/27r, constant for any S ranging 
from to 27r. Averaging over such ensemble restores rotational invar iance. In particular, for example 

{F)s = / dSF{0^5)W[5] = 1, (3.4) 
Jo 

Below, we will consider various observables, such as {N{Y^r^0^5))s as well as two-dipole correlator 
(A/'(y, ri, ^1, 5) N{Y^ r 2 -,02 -,5)) 5 averaged with respect to 5 with the weight W[S\. 

B. Single configuration solution 

In Fig. [3] we present our numerical solution of the BK equation with the initial condition (|3.2|3.3p . 

The main qualitative feature of the evolution is quick isotropisation even on a single initial configuration, without 
the ensemble averaging. One way to quantify the effect is to focus on the saturation scale, defined in a standard 
manner 

N{Y,Rs,e) = 1/2. 

The resulting saturation radius Rg is now both rapidity and angle dependent. The initial strong angular dependence 
of Rs completely disappears after evolution by about five units of rapidity (Fig. H]). 



7 




FIG. 4: Saturation radius as a function of angle and rapidity. 
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FIG. 5: a(Y) and A(Y) 



In Fig.(j5j) we plot another measure of anisotropy 

MY r) = N{Y,r,0) - NiY,r,n/2) 

^^^'"^ - NiY,r,0) + N{Y,r,7r/2) ^"^'^^ 

as a function of Y and r. Again we observe an exponentially fast disappearance of the anisotropy. 

Fig. ([5]) has a curious feature, which may or may not be important. The anisotropy A seems to be maximal at a 
fixed scale r :^ Vmax =0.5 GeV~^ independently of rapidity. The origin of this scale is not clear to us and it may just 
be a numerical accident related to the form of our initial configuration. We have checked however that A remains 
maximal at r Vmax even when R of the initial condition is varied, thus possibly hinting at another origin. At any 
rate one can follow the weakening of anisotropy by following the ratio A at rmax- We find 

A{rmax) - e-^^^ , Xa ^ 0.6. (3.6) 

We also fit in the pre-saturation regime N{Y^rmaxj^) ^ e^*^^ while N{Y^rmax^7r/2) ~ e^'^^fss]. 

We have also looked at the rate with which the scale of fixed anisotropy shrinks with rapidity. We follow the scale r 
at which A takes a constant value, say 10%. Solving for r) =0.1 leads to the scale r = a{Y). The scale r moves 
towards smaller dipole sizes with what looks like a constant rate. Our fit gives a{Y) c:^ 8.5 — Y. Note, however, that 
the fit is limited to not very high rapidities. At larger rapidities the angular dependence of N{r) is washed away and 
the anisotropy drops below 10% level for all dipole sizes. This is clearly seen on the A plot. 

It is important to stress that the wash-away of angular anisotropy occurs even for very small dipole sizes, where 
the evolution is governed entirely by the BFKL dynamics. The mechanism behind this fast isotropisation therefore 
must be rapid angular decorrelation of emitted gluons inside the BFKL ladder. 
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FIG. 6: The angle independent correlator A. 




FIG. 7: Averaged saturation radius and the fluctuation A(y, (Rs)) 



C. Averaged fluctuations 

Moving on to observables averaged over the whole ensemble of initial conditions, we first plot the fluctuation of the 
simplest angle independent correlator (Fig. [6]). 



..y . _ V{N{Y,r,0,Sr)s - {N{Y,r,e,8))l 

^ ' ^ - {N(Y,r,e,6))s ■ ^^-^^ 

We again observe the appearance of a rapidity-independent maximum at the scale Vmax^ although this maximum 
is quite shallow. The fluctuation A decays exponentially fast with rapidity with the same exponential A^^: 

A(y,w) ~ e-^-^^. (3.8) 

We have found that the exponent emerges in several other observables we looked at. Defining the angular- averaged 
saturation scale {Rs) and plotting the fluctuation A at the scale {Rs) (Fig. [7j) we again find 

A(F,{i?,)) ~ e-^-^. (3.9) 

and conclude that all fluctuations are rapidly washed away at the saturation scale. 



D. Angular Correlations 

Finally, we looked at angle dependent fluctuations of the dipole amplitude and related quantities. The first quantity 
we plot (Fig. [8]) is the correlator of two saturation scales {Rs{0i)Rs{02))5 and the normalized correlation of the 
saturation radii 



AnAY,r,9) = {R.iY,e„S))j ' 6 = 6,-02. (3.10) 

Both quantities become angle independent when evolved by about five units of rapidity. 
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The angular correlations of the dipole amplitude itself behave in a similar fashion. We plot (Fig. ^ the correlator 
{N{Y, r, 9i) N{Y, r, e2))s- and the normahzed fluctuation (Fig. 1101) [34] 



{NiY,r,0,,S)N{Y,r,e2,S))s - {NiY,r,e,,S))s {NiY,r,02,S))s 
{N{Y,r,euS))l 



61-62 (3.11) 



The normalized fluctuation decreases with rapidity approximately as 



-2\aY 



A,(F,i?,(r),^) 

Finally we plot the normalized correlation at the saturation scale as a function of rapidity (Fig. [TT]) . 
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FIG. 10: Normalized angular correlations A^. Left: y = Yo ~ 4.6; Middle: Y = 6; Right: Y = 7.5 
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FIG. 11: Angular correlations at the saturation scale Rs{Y). Left: {N{Y, Rs{Y),Oi)N{Y, Rs{Y),02)); Right: AeiY, Rs(Y),0) 



IV. POMERON LOOPS RULE. 



The question we have to ask is how do our numerical results of the previous section sit together with the picture of 
color electric fields having correlations on the scale of the saturation momentum, described in Sec. 2. Also, it is at first 
sight surprising that we do not find correlations in the target wave function which have been discussed analytically 
and numerically in [l2| (albeit for a single dipole target). 

To understand why this is the case, we note again that our calculation is done in the framework of the projectile 
dipole model, while that of [12] in the target dipole model. These two dipole models are large Nc limits of the 
JIMWLK and KLWMIJ evolutions respectively, and in the rest of this section we will not distinguish between the 
corresponding evolution and its large Nc Hmit. 

As discussed in [13] the JIMWLK and KLWMIJ evolutions are very different as far as the target wave function is 
concerned. KLWMIJ is the normal perturbative evolution of the target wave function. It is well known that in such 
evolution the number of gluons grows exponentially fast. The gluon density, which in this regime satisfies the BFKL 
equation, depends on the total rapidity as 

fl(p,F)oce^">^ (4.1) 

where c is a number of order one and Y is the total rapidity by which the target wave function has been boosted from 
'"rest". The exponential dependence is also true for the differential gluon density at any rapidity r] <Y, namely (we 
suppress the momentum dependence in the prefactor for simplicity) 

±g{p,^)c^e'^-'^e{Y-v). (4.2) 

Gluons in the wave function which are separated in rapidity by no more than 5r] ~ O(q^s) are correlated. Since the 
gluon density in the wave function exponentially grows with rapidity, any projectile that probes such a wave function 
effectively feels only the gluons in the last rapidity '"bin"'. Thus such a probe is sensitive to any correlation that 
exists between the softest gluons in the wave function. 

Now consider JIMWLK evolution. The rapidity dependence of gluon density generated by the JIMWLK evolution 
is very different. As discussed in [l3|, while the probability to emit an additional gluon in the dilute regime is 
proportional to the number of gluons in the wave function, in the dense regime this probability approaches a constant. 
The amplitude of the emission depends on the color fields in the target Ei, roughly as 

A cc ^E, ; D, = d,- gE, (4.3) 

In the dilute system this is proportional to the chromo-electric field E^ while in the dense regime, where D ^ E 
this is a constant. The evolution with constant probability of emission generates gluon density which is uniformly 
distributed in rapidity. This can be thought of as a random walk in color space as in ^] . So now one has 



(4.4) 
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with C a function of transverse momentum, but not of rapidity. It is still true, like in the KLWMIJ case, that gluons 
separated by a small rapidity interval are correlated, while the correlation disappears for gluons at very different 
rapidities. Now, as opposed to KLWMIJ however if one scatters any projectile on such a target, the projectile will 
sample gluons at all rapidities equally. Thus, for example, if the projectile consists of two dipoles, the two dipoles 
will most likely scatter on color field components (gluons) with very different rapidities. Since such fields are not 
correlated, the two dipoles will scatter independently and the two dipole scattering amplitude will not exhibit any 
correlations. 

Another way of understanding the difference between the nature of JIMWLK and KLWMIJ evolutions is the 
following. In KLWMIJ evolution one starts initially with the target wave function which contains a small number of 
gluons. The probability to emit an additional gluon in one step of the evolution is small (0(0^5)), however if a gluon 
is emitted, the end configuration strongly differs from the initial one, since the number of gluons have changed by a 
factor of order unity. Such evolution thus generates a very rough ensemble of target field (gluon) configurations. In 
this ensemble configurations with very different properties are present, albeit with small weight. Such an ensemble 
must exhibit large fiuct nations in a variety of observables. 

On the other hand, in JIMWLK evolution the probability to emit an extra gluon is large - of order unity. However 
emission of an extra gluon produces a new configuration which is hardly distinct from the existing one, since the one 
extra gluon is produced on the background of ^ gluons which already exist in the wave function. Thus the JIMWLK 
ensemble is very different - it contains many configurations, but all these configurations have very similar properties. 
The fiuctuations in observables in such an ensemble are very small. 

Thus we conclude that KLWMIJ evolution must produce (and preserve) correlations (and fiuctuations) between 
scattering amplitudes of two projectile dipoles. These are indeed the correlations discussed in [12]. On the other 
hand JIMWLK evolution must lead to disappearance of any correlations initially present in the target ensemble, as 
we have seen in our calculations. 

Given that the two approximation to high energy evolution lead to such qualitatively different answers for the 
observable we are interested in, naturally one should ask which one of them, if any should be used in quantitative 
calculations. Naively one might think, that since we are interested in high multiplicity situation, the target is dense 
and JIMWLK evolution is more suitable. This however is not the case. Recall that we expect the correlations to 
arise due to scattering at close values of impact parameter. The transverse distances in question should be smaller 
or of the order of the saturation radius Rg. This means that the scattering occurs on the components of the target 
color field with transverse momenta p < Qs- However at these momenta the target wave function by definition is still 
dilute. Even for a dense target JIMWLK evolution is not appropriate for all wavelengths. Referring to eq. ()4.3p we see 
that the amplitude is independent of the field only when we can neglect the derivative relative to the field gE in the 
expression for covariant derivative, in other words only for momenta greater than the saturation momentum. Thus 
in the dense target, we expect the gluon density to behave roughly as 

^g{p, 7?) cx e'="»''0(p - Q,(F)) + C e{Qs{Y) - p) . (4.5) 

At large momenta p > Qs the appropriate evolution is KLWMIJ and not JIMWLK. 

Thus to be able to describe fiuctuations one needs to be able to evolve the high momentum modes according 
to KLWMIJ while low momentum modes according to JIMWLK. In fact, the situation is even more complicated. 
One expects that the largest contribution to correlations comes from the modes with p of order Qs, since in the 
relative distances of order Rg should dominate the integral over impact parameter. For these modes however, the 
Pomeron merging and splitting (JIMWLK and KLWMIJ contributions) are of equal importance. In other words, this 
momentum range is sensitive at the leading order to the Pomeron loop contributions. We thus conclude that Pomeron 
loops need to be included for proper treatment of this question. 
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